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Abstract 

We propose a thermodynamic version of the Axelrod model of social influence. In one- 
dimensional lattices, the thermodynamic model becomes a Potts model of several coupled chains 
with a site (agent) interaction that increases with the site matching traits. We analytically calcu- 
late thermodynamic and critical properties for a one-dimensional system and show that an order- 
disorder phase transition only occurs at T = independent of the number of cultural traits q and 
features F of the agents. We find that the parameter q does not induce any transition or anomaly 
in the thermodynamic model, as it does in the standard social model that violates detailed balance. 
The one-dimensional thermodynamic Axelrod model belongs to the same universality class of the 
Ising and Potts models notwithstanding the increase of the internal dimension of the local degree 
of freedom (agent). 

PACS numbers: 89.20.-a, 89.75.Hc 
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I. INTRODUCTION 



Statistical mechanics has helped us to understand macroscopic collective phenomena of 
a system in terms of its microscopic constituents and their interactions. Sociophysics is 
born under the idea of transferring this premise to social phenomena (formation of public 
opinions, actions of social groups, epidemic spreads, etc.), motivated by the fact that in 
many cases those phenomena observed in complex networks are analogous to those found in 
systems in thermal equilibrium [1]. 

The Axelrod model is an example where a local (microscopic) rule is used for the study of 
global social behaviors [2]. This model studies, in particular, the dissemination of cultures 
among interacting individuals or agents in which (a) the more culturally similar the agents, 
the greater the chance of interaction between them, and (b) interaction increases similarity 
between individuals. More explicitly, the model considers that an agent located at the i th 
site of a lattice is defined by a set of F cultural features (e.g., religion, sports, politics, etc.) 
represented by a vector Oi = (an, a i2 , ■ a^p). Each feature a^ can take integer values in 
the interval [l,q], where q defines the cultural traits allowed per feature and measures the 
cultural variability in the system. There are q F possible cultural states. 

The model's dynamics is as follows: (1) Choose randomly two nearest neighbor agents % 
and j, then (2) calculate the number of shared features between the agents £y = J2k ( W>ffjk- 
If < £ij < F, then (3) pick up randomly a feature k such that ^ a^ and with 
probability £ij/F set a^ = ajk- These time steps are iterated and the dynamics stops when 
a frozen state is reached; i.e., either £y = or = F, Vz,j. A cluster is a set of connected 
agents with the same state. Monocultural or ordered phases are composed of a cluster of 
the size of the system where £ij = F,Wi,j. Multicultural or disordered phases consist of two 
or more clusters. 

One of the main features of this model is a change of behavior, at a value q c , from a 
monocultural state, where all agents share the same cultural features, to a multicultural 
state, where individuals mostly have their own features, as the cultural diversity q increases 
j^l . This change can be characterized by an order parameter <p that is usually defined as the 
average size of the largest cultural cluster C max normalized by the total number of agents 
N in the system; = C max /iV. In the monocultural (ordered) state — > 1 and in the 
multicultural (disordered) state — > 0. 
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This change of behavior in the Axelrod model is often associated to phase transitions 
that occur in some physical systems 3|-|5|. However, it is known that true phase transitions 
occur only in the thermodynamic limit. The suggested "transitions" when noise is added 
to the Axelrod model, depend on finite size effects and the critical values q c disappear in 
the thermodynamic limit as q c cx 1/N (5), 1^|. On the other hand, social systems are 
believed to be out of equilibrium. The association of the phenomena in the Axelrod model 
with physical systems is thus purely figurative. Nevertheless, concepts like temperature, 
critical phase transition, critical point, applied magnetic field, among others, have previously 
been used in studies of the Axelrod model (5), 8]. Here, with the aim of finding a possible 
thermodynamic role of the parameters F, q, and q c , we introduce a thermodynamic model 
that captures the main characteristics of the nonequilibrium Axelrod model. 

For our model, hereafter called thermodynamic Axelrod, we use the interaction rule or 
number of shared features of the Axelrod model to construct a new Hamiltonian dis- 
tinguishable from the F parallel one-dimensional (ID) Potts models in that the interaction 
strength between agents increases with £ij. This feature of the interaction precipitates or- 
dering preempting fluctuations. In the thermodynamic Axelrod model F is related to the 
coupling energy of the system and q has the same meaning as in the Potts model. 

The paper is structured as follows: In Sect. II we demonstrate that the standard Axelrod 
model does not satisfy the detailed balance condition. In Sect. Ill we present the mathe- 
matical framework for computing the physical properties of our model, which are actually 
determined in Sect. IV. We compare the original and thermodynamic Axelrod models in 
Sect. V and give some conclusions in Sect. VI. 



II. AXELROD MODEL: OUT OF EQUILIBRIUM 

Before we introduce the thermodynamic version of the Axelrod model, here we demon- 
strate that the original version does not satisfy equilibrium conditions by showing that it 
does not obey detailed balance. 

Let a link between two sites i and j be of type n if they share n components = n), 
P n be the probability that the system is in a state with links of type n, and W nm be the 
transition probability per unit time from a state with type-n to one with type-m links. W nm 
being time independent. Since in the dynamics of the Axelrod model a feature k is changed 
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with probability d-ij/F to make two sites have one more component in common (crjfc = cr^), 
we have 



. % for m = n + 1 
W nm = { F . (1) 

otherwise 



The master equation 



alP 

T = H l W nmPm(t) ~ W mn P n (t)} 



dt 

controls the time evolution and determines the equilibrium states of the system 9j. As is 
well known, an equilibrium state is characterized by a probability P n independent of time t. 
Thus, 



n . 



Wn m P m - W mn P n ) = for any 

rn 

The detailed balance relation implies that 

n,m. (2) 

In the Axelrod model Eq. (j2J) cannot be satisfied since one side or the other is always zero 
according to Eq. (JT]). This feature of the model introduces some very strong constraints into 
both the evolution and the time- independent states of human interactions that emphasize 
nonequilibrium; i.e. i) completely different individuals do not interact, ii) individuals "con- 
form" once they have modified their cultural profile, and iii) individuals that are alike, that 
interact, increase their similarity at interaction. These rules yield frozen final states with no 
further evolution or fluctuation. 

The thermodynamic model we propose relaxes all of the previous constraints, while pre- 
serving similarity increasing interactions and introducing fluctuations controlled by a tem- 
perature parameter. These conditions permit arriving at a dynamical equilibrium state in 
the regular sense of statistical mechanics. On the other hand, we are interested in evaluating 
whether behaviors reported in nonequilibrium networks models survive a thermodynamic 
driven scenario. 
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III. NATURE OF THERMAL SOCIETY OF AGENTS 



Two basic questions come to mind when proposing a thermodynamic model for society, 
namely: what is the meaning of temperature and detailed balance? In the Axelrod model 
each agent is described as a vector whose F components can take q values that reflect 
the variety of postures an individual can have on a particular scope of action in society. 
Temperature may be thought as related to an energy scale that competes with the regular 
interactions between individuals; a high temperature renders their interactions moot, while 
a low temperature leads to a domination of the individual interactions and to the system 
settling into a minimum energy state. In this sense, in a society at high temperatures 
individuals would have fluctuating positions with a small relation to those of their circle 
of interaction, whereas at low temperatures individuals would pay much attention to their 
circle of interaction and tend to lower posture differences. In terms of a statistical ensemble 
approach, a group of individuals of the society would be subjected to a temperature reservoir 
that competes with individual interactions. 

What sets this temperature scale might be any motive for confusion, for speculation or 
for uncertainty that could disrupt bonds between individuals. This would be consistent with 
regarding temperature as a parameter coupled to entropy, in the sense of the First Law of 
Thermodynamics. In social terms, as the temperature increases the agents communicate less 
effectively, their coupling decreases, they are less convincing in a sense to their neighbours 
who choose or preserve their own positions without regard to their peers. This has entropic 
benefits since there is an increasing amount of ways people can disagree. On the other 
hand, when the system approaches to the critical point each individual takes the states of 
its neighborhood and each neighborhood takes the state of its around neighborhood and in 
this way all individuals remain correlated with the rest of individuals of the system. 

Thermal equilibrium is a situation where some Helmholtz-type free energy is minimal, 
and energy fluctuations exist subject to the condition of detailed balance. From Eq. (J2J) 
one can see that detailed balance dictates that an individual in a highly probable state n 
should balance with an individual in a less probable state m. It is reasonable to expect 
that it is easier for the individual in state m to adjust itself to the mainstream rather than 
the other way around. This can be seen as peer pressure or pressure to conform to the 
norm. This is an interesting perspective in the sense that varying temperature can yield 
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thresholds for phase changes and set tipping points for collective behavior. On the other 
hand, a single temperature may not be set for different scopes of action, since a religious 
posture is certainly less fluctuating than a political posture. In our model this may be 
taken into account by weighting the interaction factor J^- (see below) that depends on the 
cultural-featured F-vector. We will not address this last point in the present work. 

IV. THERMODYNAMIC AXELROD MODEL 

To mimic the interaction rule of the Axelrod model, in which the interaction probability 
is proportional to the number of shared features, the Hamiltonian for the thermodynamic 
model is defined as 

F N ( H \ 

k=l ij ^ ' 

with the interaction factor 

F 

Jij = y]J5((Ti n , a jn ) . (4) 

n=l 

Here, H k works as an applied magnetic field that can point in one of the Potts-model-like 
directions k, can take values 1, q, and has an energy weight proportional to the magnitude 
H. (jjfe = 1, . . . , q specifies each of the F variables (an, a i2 , ■ ■ ■ , cr iF ) of the agent at at the ith 
lattice site. iV is the size of the system. The second term in the Hamiltonian is symmetrized 
for convenience. 

This Hamiltonian is evidently inspired on the Potts model, with the significant distinction 
that the interaction factor always depends on the global state of the F-vector and not 
on the state of the particular Potts variable. This is a somewhat rare type of Hamiltonian 
interaction, in a sense similar to a non-linear sigma model where a vector interaction occurs 
subject to no—ion of the interacting vectors Q. Our .node, on a !D .attics is like au 
F-layer coupled Potts model; it is a quasi one-dimensional system so that some signatures 
of the two-dimensional Potts model are expected as crossovers. 

We consider a ID chain of iV sites occupied by g F -valued agents and use the transfer 
matrix method to compute the relevant physical properties. For periodic boundary condi- 
tions cx(jv +1 )fc = oifc, the partition function corresponding to the above Hamiltonian can be 
expressed as 
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z = EE-EII«p 

C\ (72 (TJV i=l 



^ ( /3Ji(i+i)5((Tik, a 
k=i ^ 

<5((7jfc, iffc) + 5((T(j + i) fc , iffc 



(*+!)&; 
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Here, (5 = l/fc^T and we introduced & q x q transfer matrix W with elements 



(5) 



(ai\W\a i+1 ) = exp 



F ( 

fe=l ^ 



(6) 



The eigenvalues Xj of the transfer matrix are determined from the solution of the secular 
equation DetlW 7 — XE\ = 0. Then, the partition function can be written as 



Z — Aj v + A J 2 V + • • • \q F — X ! ma x I 1 + — I 1" T/f" I • 

where X max is the largest eigenvalue. In the thermodynamic limit (N — > oo), Z = X 



(7) 



A' 



Then, a standard procedure 
properties. 



Ill ] can be followed to obtain thermodynamic and critical 



V. CASE F = 2, q = 2 

Here, we present explicitly the simplest non-trivial case F = 2, q = 2, as an explicit 
illustration. However the analytical calculations were performed for higher values of F and 
q. The transfer matrix takes the form 

/ e 2/3^H+il3J e lP^H+pj e lfDnH+pj e /3pH \ 

w = l (8) 

The largest eigenvalue for this matrix is 
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where 



A, 



D A/3J+l3tiH 



+ e 



4I3J+2/3/1H 



+ e 4/3J + e 



3 [b + V4a 3 + b 2 \ 



1/3 



1 



3 2V3 L 



b + y/Aa? + 6 2 



1/3 



(9) 



fl = _g e 2/3J+/3 Mj H- _|_ ^J+PfiH _|_ e 8/3J+/3 M H _ 2e 4/3J+2/3 Mj H' _ g e 2/3J+3/3 M # 
j^ e Al3J+WvH _|_ e Spj+3/3nH _ e 8/3J+4/3fj,H _ g 8/3J _ ^PfiH 



+lge 6 / 8J+2 / 3/iH _ g e 8/3J+2/3 M H _ g e 12/3J+2/3 Mj H- + ^Qg^J+S^/xH _ ^4/^+3/3^ 
_72 e 6/3J+3/?M^ _|_ 18e 8/3J+3/3^ + -^12/37+3/3/^ + lge 2/3 + g g 4/3 

+18e 6/3J+4/3AlH - 6e 8/3J+4/3AlH - 6e 12/3J+4/3 ^ + i8e 6/3J+5/3AlH - 2>e^ J+hPliH 

_ 3e 12/3J+5/3^ _|_ 2e 12/3J+6/3^ + 2 e 12/3J - 16e 3 ^ H . 

The free energy, magnetization, magnetic susceptibility, and specific heat are given in 
terms of X max : 

F=-Nk B T\n\ max , (10) 



dT Nk B Td\ max 

M = = — , (11) 

dH A dH 7 



= m_ d ( Nk B Td\ max 

X air azj 1 \ air J ' v / 



dH dH \ \ max dH 



^d 2 ? Nk B Td\ 

max , Ar ; ml 

d ( i d\ 

max \ / 1 o \ 

c = ~ T df^ = 2 ^^r + t [x^^r) ■ (13) 

Figures dH show 3D plots of the magnetization, magnetic susceptibility, and specific heat 
in terms of k B T/J and fiH/k B T. All these properties behave qualitatively as they do in the 
Ising and Potts models in ID lattices, no new qualitatively striking behavior is obvious. 
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FIG. 3: Specific heat of the ID thermodynamic Ax- 
ehod model for the case F = 2 and q = 2 

The magnetization (Fig. [I]) goes to as H — >• at any finite T. At T = the magne- 
tization saturates to its maximum value for any H. This implies a spontaneous transition 
to an ordered state only at T = 0. As T becomes finite, the magnetization saturates to its 
maximum only at large H. The susceptibility (Fig. [2]) diverges as T — > and H — > 0. For 
zero field the specific heat goes through a maximum before vanishing smoothly as T = 
(Fig. [3]). Thus, the differences in internal symmetry and interaction in the thermodynamic 
Axelrod model with respect to the Ising model do not seem to be relevant. We go back this 
point later in the discussion on the critical exponents. 
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In terms of the thermal society of agents, the above results imply that the fluctuating 
positions will be always so large as for the individuals to conform spontaneously to a par- 
ticular posture, unless there exist a highly persuasive agent or mass media. On the other 
hand, the specific heat can be regarded as the ability of the society to reject the increase of 
spontaneous changes of posture. This ability decreases as the occurrence of these changes 
goes down and disappears for zero occurrence. 



A. Spatial correlations 



We now calculate the correlation function 



G(i,l+j) = ((Ti<Ti+j) - (<7i)(<7;j) 



(14) 



using the transfer matrix method [121 ] . The two-agent correlation function is given by 

(ww) = ^tr[AW j AW N ~ j ] , 



(15) 



.4 



(16) 



where A = \o'i)a i (o'i\. In the original space for F = 2 and q — 2, the state matrix 

( I o\ 
2 
2 
\0 4/ 

Following a standard procedure, we evaluate Eq. (1151) in a basis where W is diagonal. For 
H = 0, the unitary matrix that diagonalizes W is 

/-l 1 l\ 
-1-11 
1-11 

V i o 11/ 



(17) 
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The evaluation yields 
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Here, A4 and Ai are the largest and the second largest eigenvalues, respectively, and the 
evaluation was performed in the thermodynamic limit. 

The average of the agent at site i, (<7j), is evaluated in the same manner: 



(a t ) = -tr[AW N ] 



9 
4 



Then, 



(19) 



(20) 



The correlation function has the same form as in the Ising and Potts models 12j, ll3], with 
the correlation length 

(21) 



Since Ai 



ln(A 4 /A 1 )' 

-1 + e 4J / fcsT and A 4 = 1 + 2e J//fesT + e 4J//fcsT , the correlation length becomes 



1 



(22) 



In [(1 + 2e J / k B T + e 4J / fc s T ) / (-1 + e 4J / fe ^)] " 

As in Ising and Potts models, for H = the two largest eigenvalues become degenerate at 
T = 0, leading to a divergence of the correlation length and, therefore, to a zero-temperature 
phase transition. 

In social terms, the correlation length measures the distance at which there are relations 
between agents beyond their own mean values. This is a causal or influence relationship 
in the sense that changing the opinions in one place generates an influence that causes 
change up to the correlation length. This influence operates through local interactions. In 
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general, correlation lengths are linked to a return force that limits their growth. In terms 
of social influence the existence of a correlation length invokes a limit to the propagation of 
influence, some return force avoids the propagation of fluctuations beyond a certain distance. 
When no return force is present fluctuations run wild and influence runs over the whole 
society. The return force is a measure of the cost of producing the fluctuation. If there is no 
cost, fluctuations diverge and a change of opinion in one place generates a fluctuation that 
propagates throughout the system. 



B. Critical exponents 

To get the critical behavior of the thermodynamic properties one needs to evaluate them 
near the transition temperature T c . Usually, for this purpose T is replaced by the parameter 
t = (T — T c ) /Tr . However, for ID models the transition occurs at T c = with exponential 
singularities In this case, the regular parameter t as defined is inappropriate. A different 



critical point approach parameter t = e _A//fcsT ll| is required in the ID models to convert 
the exponential singularities in T into power-law singularities in t. The constant A has been 
thus far taken arbitrarily, which leads to ambiguous critical exponents. 

Here, we propose that A be given by half the energy difference between the ground 
state and the first excited state of the system. In this way A is no longer arbitrary. This 
convention correctly unifies, independent of the interaction strength, the Ising and Potts 
exponents. For the ID Ising model the energy difference between the ground state (say, all 
spins aligned up) and the first excited state (one spin aligned opposite to the others) is 4J; 
then, for this case A = 2J. This agrees with the choice of p = 2 in t = e~ pJ//fcsT to bring 
together the exponents of the ID discrete-symmetry models. 

For the Potts and thermodynamic Axelrod models, the energy of the ground state is 
—JNF 2 and of the first excited state is -J [(N - 2)F 2 + 2(F - l) 2 ]. Then, 

A = J(2F - 1) . (23) 

In the case of the Potts model, F — 1 and A = J. This value of A yields critical exponents 
of the Potts model that agree with those of the Ising model For the thermodynamic 
Axelrod model we are analyzing, F = 2 and A = 3J. 

We are now prepared to estimate the critical exponents of our model for the case F = 2 



12 



and q — 2. We define h = fiH/k B T and have t = e ~ 3J / fcsT . For H = and t — > the 
singular part of the free energy, Eq. ( ITUl) . for a zero temperature transition [ill] becomes 



/(i) = £±i^ t . (24 ) 



Since / ~ t a , a — 1. In the same limits, the magnetization, Eq. (fllj) . is 

This means from m ~ t p that /3 = 0. Now, for t = and -)• 0, the magnetization becomes 

m(h) = ^ ~ 1 . (26) 
Since m ~ /i 1 /" 5 , the exponent <5 — >■ oo. The low- field susceptibility is obtained from Eq. ( fl2l) 

«.M = ^~«-'. <") 

The susceptibility should go as t -7 , then 7 = 1. The specific heat, Eq. ( flBl) . for = and 
t is 

Then, from C ~ t Q one gets a = 1. Finally, from Eq. f l22l) the correlation lenght 

£(f)~f -1 . (29) 

The correlation length goes as t~ u ; then, u = 1. 

In conclusion, the differences in the internal space dimensionality (F = 2) and interaction 
strengths do not alter the Ising universality class of our model. 



VI. THERMODYNAMIC AND SOCIAL AXELROD MODELS: A COMPARISON 

The original Axelrod model in the absence of external field displays a crossover from a 
disordered to an ordered state at certain value q c that depends on dimensionality 
Moreover, the way this crossover takes place depends on the cultural feature F; for F > 2 
it is abrupt and for F = 2 it is smooth 3]. For F = 2 the crossover occurs for q < 5 [4j. As 
shown above, our thermodynamic Axelrod model in ID has a continuous phase transition to 
an ordered state at T = 0, just as expected for discrete symmetry systems and as is found 
in the Ising and Potts models. To see if there is any correspondence between the trait q and 
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feature F variables in the two Axelrod models; that is, if q and F have any effect on the 
thermodynamic model that can be mapped into the social model, we carried out analytical 
calculations for different F and q, sweeping the latter across the region were the crossover 
happens in the social model. 

Figure @] depicts the analytical results for the temperature dependence of the order pa- 
rameter for F = 2 and q = 2,3,4,5. For comparison purpose only, we also show curves 
for F = 3 and q — 2, 3. A monotonic and smooth dependence on q is observed for F = 2, 
indicating no signal of crossover or phase transition at all. For both values of F, the order 
parameter behaviors are similar to those seen in the Potts model [3] . We note that for the 
analyzed values of F as q increases the transition becomes gradually sharper, which could 
mistakenly suggest a change in the transition type as a function of q. We know that for ID 
systems this is not allowed and that the observed effect is due to numerical limitations as 
H -»• 0. The data of Fig. H were obtained for H = 0.005. 

The results suggest that thermodynamically the ID Axelrod hamiltonian does not exhibit 
a phase transition as a function of the variable q. The transition depends only on tempera- 
ture and in one-dimension it takes place in T=0. It is interesting that in the non-equlibrium 

n n 

model with noise according to Klem et. al [14J, [15(, if one equates temperature to noise, 
for any finite noise intensity the system is "polarized" for large systems, in the sense that 
it breaks up into finite islands. This is equivalent to our zero magentization phase for any 
finite temperature. 



VII. SUMMARY 



We presented a thermodynamic counterpart of the Axelrod model of social influence. For 
the thermodynamic version we calculated thermodynamic and critical properties in ID and 
showed that an order-disorder phase transition occurs at T = 0, independent of the cultural 
trait q and feature F variables. A relevant result is that for our thermodynamic model the 
parameter q does not induce any transition or anomaly in the system, as it does in the social 
model. We showed that the thermodynamic Axelrod model belongs to the same universality 
class of the Ising and Potts models. 
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FIG. 4: Magnetization or order parameter of the thermodynamic Axelrod model. The group of 
curves to the left correspond to F = 2 and the one to the right to F = 3. No anomaly is observed 
as q varies from 2 to 5 for the F = 2 case. In the original Axelrod model a crossover from an 
ordered to a disordered state is seen below q = 5 in the one-dimensional case [J]. 
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